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THE EFFECTS OF THERMAL CONDUCTION ON THE HOT ACCRETION FLOWS 
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ABSTRACT 

Thermal conduction has been suggested as a possible mechanism by which sufficient extra heating is 
provided in radiation-dominated accretion flows. We consider the extreme case in which the generated 
energy due to the viscosity and the energy transported by a saturated form of thermal conduction are 
balanced by the advection cooling. For the steady-state structure of such accretion flows a set of self- 
similar solution are presented. Based on these solutions while the radial and the rotational velocities are 
both sub-Keplerian, increasing the level of thermal conduction has the effects of decreasing the rotational 
velocity, but increasing the radial velocity. Conduction provides extra heating and the temperature of 
the gas increases with thermal conduction. 

Subject headings: accretion, accretion disks - hydrodynamics 



1. INTRODUCTION 

The thin accretion disk model describes flows in which 
the viscous heating of the gas radiates out of the system 
immediately after generation (Shakura & Sunyaev 1973). 
However, another kind of accretion has been studied dur- 
ing recent years where radiative energy losses are small so 
that most of the energy is advected with the gas. These 
Adection-Dominated Accretion Flows (ADAF) occur in 
two regimes depending on the mass accretion rate and the 
optical depth. At very high mass accretion rate, the op- 
tical depth becomes very high and the radiation can be 
trapped in the gas. This type of accretion which is known 
under the name 'slim accretion disk' has been studied in 
detail by Abramowicz et al. (1988). But when the accre- 
tion rate is very small and the optical depth is very low, we 
may have another type of accretion (Narayan & Yi 1994; 
Abramowitz et al. 1995; Chen 1995). However, numeri- 
cal simulations of radiatively inefficient accretion flows re- 
vealed that low viscosity flows are convectively unstable 
and convection strongly influences the global properties 
of the flow (e.g., Igumenshchev, Abramowicz, & Narayan 
2000). Thus, another type of accretion flows has been 
proposed, in which convection plays as a dominant mech- 
anism of transporting angular momentum and the local 
released viscous energy (e.g., Narayan, Igumenshchev, & 
Abramowicz 2000). 

This diversity of models tells us that modeling the hot 
accretion flows is a challenging and controversial problem. 
We think, one of the largely neglected physical ingredi- 
ent in this field, is thermal conduction. But a few authors 
tried to study the role of "turbulent" heat transport in 
ADAF-like flows (Honma 1996; Manmoto et al. 2000) 
Since thermal conduction acts to oppose the formation 
of the temperature gradient that causes it, one might ex- 
pect that the temperature and density profiles for accre- 
tion flows in which thermal conduction plays a significant 
role to appear different compared to those flows in which 
thermal conduction is less effective. Just recently, Menou 
(2005) studied the effect of saturated thermal conduction 
on optically thin ADAFs using an extension of self-similar 
solutions of Narayan & Yi (2004). His solutions suggest 



that thermal conduction may be an important physical fac- 
tor to understand hot accretion onto dim accreting black 
holes. 

The aim of the work presented here is to investigate 
whether thermal conduction can affect the dynamics of 
an optically thick radiation-dominated accretion flow. In 
§2, we present a height-integrated set of equations of the 
model. Self-similar solutions are investigated in §3. The 
paper concludes with a summary of the results in §4. 

2. FORMULATION 

In order to implement thermal conductivity correctly it 
is essential to know whether the mean free path is less 
than (or comparable to) the scale length of the temper- 
ature gradient. For electron mean free path which are 
greater than the scale length of the temperature gradient 
the thermal conductivity is said to 'saturate' and the heat 
flux approaches a limiting value (Cowie & McKee 1977). 
But when the mean free paths are much less than the 
temperature gradient the heat flux depends on the coeffi- 
cient of thermal conductivity and the temperature gradi- 
ent. Generally, thermal conduction transfers heat so as to 
oppose the temperature gradient which causes the transfer. 
Menou (2005) discussed hot accretion likely proceed un- 
der weakly-collisional conditions in these systems. Thus, 
a saturated form of "microscopic" thermal conduction is 
physically well-motivated, as we apply in this study. How- 
ever, one of the primary problems for studying the effects 
of thermal conduction in plasmas is the unknown value of 
the thermal conductivity. 

We consider an accretion disk that is axisymmetric and 
geometrically thin, i.e. h/r < 1. Here r and h are, respec- 
tively, the disk radius and the half-thickness. Our model 
generalizes the usual slim disks around a neutron star or 
a black hole (e.g., Muchotrzeb & Paczyhski 1982; Mat- 
sumoto et al. 1984; Abramowicz et al. 1988) by including 
the effect of thermal conduction. The disk is supposed to 
be turbulent and possesses an effective turbulent viscos- 
ity. We assume the generated energy due to viscous dis- 
sipation and the heat conducted into the volume is con- 
cerned are balanced by the advection cooling. Consider 
stationary height-integrated equations described an accre- 
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tion flow onto a central object of mass M*. In absence of 
mass outflows, the continuity equation reads 

M = -27rrw r S, (1) 

where M is the accretion rate, v r is the accretion veloc- 
ity (and v r < 0) and X = 2hp is the surface density at a 
cylindrical radius. Also, p is the midplanc density of the 
disk. 

The equation of motion in the radial direction is 



ClVr- 



dr r dr r^ 

where P = 2hp is the integrated disk pressure written, 
in a form compatible with radiation-dominated flow, as 
P ~ 2/iaT 4 /3, where a is black body constant and T de- 
notes the midplane temperature of the disk. Similarly, in- 
tegration over z of the azimuthal equation of motion gives 

pB *s <"*> - < 3 > 

where v is a kinematic viscosity coefficient and we assume 

v = ac s h, (4) 

where a is a constant less than unity (Shakura & Sun- 
yaev 1973). Moreover, h = (Cg/i^r, where c s and vk 
are sound speed and the Keplerian velocity, respectively: 
c s = v/P/S and w K = y/GM*/r. 

Now, we can write the energy equation considering the 
cooling and the heating processes of the accretion flow. 
The locally released energy due to viscous dissipation, 
Qvis, and the energy transported by thermal conduction, 
Qcond, are balanced by the advection cooling, Q a d v - Thus, 

Q&dv — Qcond ~t~ Qvis- (5) 

The advection cooling reads (Abramowicz et al. 1998) 

Qadv = ^ — j _ "F 7T71 + 4 ~ 3 ^ Ti — ' 6 

2itr A p r 3 — 1 dlnr dlnr 

where /? is the ratio of the gas to the total pressure and 

(4 -3/?) (7-1 



dP GM* 

: 



(2) 



r 3 = i 



(7) 



/?+12( 7 -l)(/3-l)' 

Since we are considering the extreme case of radiation- 
dominated accretion flow, we have (3 = and so T 3 =2/3 
Also, we have 



Q yiB = uVr 2 [±(^)] 2 , 
dr r 



and (Cowie & McKee 1977) 



2h d 2 

^cond <7~i~ V s ) 



(8) 



(9) 



where as we have already mentioned F s = 5$ s pCg is the 
saturated conduction flux on the direction of the temper- 
ature gradient (<J> S < 1). 

3. SELF-SIMILAR SOLUTIONS 

The self-similar solution is not able to describe the 
global behaviour of the accretion flow, because no bound- 
ary condition has been taken into account. However, as 
long as we are not interested in the the behaviour of the 
flow at the boundaries, such solution describes correctly 
the true solution asymptotically at large radii. We assume 
that each physical quantity can be expressed as a power 



law of the radial distance, i.e. r v , where the power index v 
is determined for each physical quantity self-consistently. 
The solutions are 



E(r)=aE (-)- 



-1/2 



v T (r) = 



GM* f r ,_ 1/2 
ro V ' 



GAL 



ro r 



(fr 1/2 > 



P(r) ^ d ^Ml ( ^ ) -3/ 2 

ro r 

h{R) = fro(-), 
ro 



(10) 

(11) 
(12) 
(13) 
(14) 



where So and ro provide convenient units with which the 
equations can be written in the non-dimensional form. 
Thus, we obtain the following system of dimensionless 
equations, to be solved for a, 6, c, d and /: 

(15) 
(16) 



(17) 
(18) 

(19) 




af -2d = 0, 




fb 2 a + 5$ s d\ 



where m = M/(27rE \/GM»ro) is nondimensional mass 
accretion rate. After some algebraic manipulations we find 
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b=\l 



a 2 / 4 



16 



(20) 

(21) 
(22) 

(23) 



and / is obtained from a forth order algebraic equation: 

160$, 



9a 2 / 4 + 20 f 2 



9a 



-f - 16 = 0. 



(24) 



We can solve this algebraic equation numerically and 
clearly only the real root which corresponds to positive 
b 2 is a physical solution. 

If we neglect the thermal conduction (i.e., $ s = 0), the 
above equation reduces to a second order equation which 
can be solved analytically (Wang & Zhou 1999; Shadmehri 
& Khajenabi 2005). In this case, the typical behaviour of 
the solutions can be summarized as follows (Shadmehri & 
Khajenabi 2005): (1) The surface density increases with 
the accretion rate, and decreases with viscosity coefficient 
a; (2) But the radial velocity is directly proportional to a; 
(3) The gas rotates with sub-Keplerian angular velocity, 
more or less independent of the coefficient a; and (4) the 
opening angle of the disk is fixed, independent of a and rh. 
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Note that Shadmehri & Khajenabi (2005) applied a diffu- 
sive prescription for the viscous stress tensor in their anal- 
ysis of the self-similar solution of optically thick advection- 
dominated flows, but Wang & Zhou (1999) applied the ap 
prescription. However, the typical behaviours of the solu- 
tions are the same, irrespective of some differences in the 
coefficients. 

The equations (20), (21), (22), (23) together with equa- 
tion (24) describes self-similar behaviour of the optically 
thick advection dominated accretion disk with saturated 
thermal conduction. The algebraic equation (24) shows 
that the variable / which determines the opening angle of 
the disk depends only on the a and § s . This behaviour is 
similar to the case without thermal conduction. Thus, one 
can deduce that the surface density and the pressure are 
directly proportional to the mass accretion rate according 
to the equations (20) and (23). But the rotational and the 
radial velocities are both independent of the mass accre- 
tion rate (see equations (21) and (22)). Also, we can say 
that the rotational velocity is sub-Keplerian according to 
the equation (21). 

The dependence of the properties of the rotational and 
the radial velocities on a and $ s is illustrated in Figure 
1. We consider three values of the viscosity parameter, 
a = 0.2, 0.02 and 0.002 for which the velocities are plotted 
as a function of $ s . The solid and dashed lines represent 
the ratios v v /vk and v t /vk, respectively. While the so- 
lution without thermal conduction are recovered at small 
$ s values, we can see significant deviation as $ s increases. 
Clearly, the rotational velocity decreases with the magni- 
tude of conduction parameter. But for the radial velocity 
the behaviour is different, i.e. the velocity increases as 
$ s increases. For a fixed viscosity parameter a, the solu- 
tion reaches the non-rotating limit at a specific value of 
$ s which we denotes by If we extend our solution 
to <!> s > the equation (21) gives negative b 2 which is 
clearly unacceptable. The critical value $g can be calcu- 
lated easily: 

^ = -J^VV4^+1~. (25) 

The above equation shows the critical magnitude of con- 
duction parameter for which the solution tends to the non- 
rotating limit depends only on the viscosity coefficient so 
that as a increases, the breakdown of the solutions occurs 
at larger values of <& s . But as Menou (2005) mentioned in 
his analysis of optically thin ADAFs with saturated con- 



duction, we think this behaviour is a simple pathological 
feature of ID height-integrated equations we solved. 

Solutions of equation (24) show that for a fixed a, the 
variable / increases with conduction parameter $ s . Since 
the gas temperature varies directly in proportion to /, we 
can say as the level of thermal conduction is increased the 
temperature of flow in comparison to solutions without 
thermal conduction increases. In fact, thermal conduc- 
tion behaves as an extra heating source. Also, equation 
(20) shows the surface density decreases as $ s increases. 
We can compare the relative values of advection cooling, 
viscous heating and the transported energy by conduc- 
tion using our self-similar solutions. One can simply show 
that Q v i s /Qadv = 2(fe//) 2 and Q co „d/Qadv = 20$ B /(9a/). 
While both these ratios don't depend on the mass accre- 
tion rate, the viscosity parameter a and the thermal con- 
duction parameter $ s determine the relative importance 
of the various physical processes in the system. Viscous 
energy dissipation decreases as $ s increases, but the con- 
tribution of the energy transported by the conduction be- 
comes more significant. 

4. CONCLUSIONS 

In this paper we have studied an optically thick 
ADAFs with thermal conduction. Considering the weakly- 
collisional nature of hot accretion flows, a saturated form 
of thermal conduction has been used as a possible mecha- 
nism of transporting of energy. Although our self-similar 
solutions have their own limitations, they can illuminate 
some possible effects of the saturated thermal conduction 
on radiation-dominated accretion flows. We can summa- 
rize the main features of our solutions: (1) The surface 
density is in proportion to the mass accretion rate and 
it increases as the level of the thermal conduction is in- 
creased. (2) The temperature of the gas increases with 
thermal conduction parameter and the conduction term 
in the energy equation acts as an extra source of heating. 
(3) The angular velocity is sub-Keplerian, but it decreases 
with increasing thermal conduction parameter so that it 
tends to a non-rotating limit. (4) However, the radial ve- 
locity increases with thermal conduction parameter. Con- 
sidering these significant effects of the thermal conduction, 
global transonic solutions representing optically thick (or 
thin) ADAFs with thermal conduction would be an inter- 
esting topic for future studies and our results are a step 
toward this goal. 
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Fig. 1. — Profile of the rotation velocity (v^/vyi, solid lines) and the radial velocity (v r /vn, dashed lines) as a function of the saturation 
constant, $ s . In order to make easier comparison the rotational and the radial velocities are scaled by the Keplerian velocity, «k- Each curve 
is labeled by corresponding viscosity coefficient, a. 



